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INTRODUCTION 


A prototype interferometric spectrometer system is being built by NASA to explore 
the potential of the technique for applications involving the visible and near- 
ultraviolet portions of the electromagnetic spectrum. As described in the paper by 
Parsons, Gerlach, and Whitehurst (1982), which will hereafter be referred to as PGW, the 
system is limited only by the frequency bandpass of the optical components used in the 
system, the quality of the optical components, and ultimately by the memory capacity of 
the computer. In PGW, however, the tradeoff between the wavenumber resolution of the 
produced spectrum, the bandpass limits of the optics, and the number of samples obtained 
from the interferogram were not delineated explicitly. The prototype Ultraviolet Inter- 
ferometric Spectrometer (UVIS) instrument is expected to be configured several different 
ways to ascertain its suitability for various applications. To exploit its inherent 
flexibility, a reference document containing the parameter tradeoffs involved would be 
very useful. This paper is designed to be that reference. 


BACKGROUND 


Ihe Interferogram 


In PGW, the transform pair linking the wavenumber, o, and path difference, x, 
domains was derived and found to be 


y(x) 



P(o)e2TriaX dg 
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dx 


(la) 
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where y(x) Is the Interferogram, or interference pattern, produced by the spectrometer 
when exposed to an Incoming light beam having the spectral shape P(cr). In actuality, 
of course, it is not possible to perform the integral calculations because the spectrum 
is too complex to permit closed-form solutions. Therefore, the integral transforms must 
in fact be replaced by a discrete Fourier transform pair. In this section, this relation- 
ship will be derived so that the pertinent parameters for interferometric spectroscopy 
and their interconnections can be described mathematically. 

Figures 1 and 2 show graphically the steps in the path difference and wavenumber 
domains, respectively, that carry the interferogram and incoming spectrum at the top to 
the recovered discrete spectrum at the bottom. Each step will be discussed in detail and 
sequentially. The starting point is the production of the interferogram y(x) by the 
optical interference produced in the Michelson interferometer. Therefore, the linkage 
between Figures la and 2a is the optics in the instrument. 

Once the path difference signal is produced, it must be captured in some way and 
this is done using an analog-to-digital converter. In the next section, the technique 
used to produce samples evenly spaced in path difference will be described in more detail. 
The sampling process can be described by the equation 

00 

y(x)a(x) = y(x) ^ 6(x-kAx) (2) 

k=-eo 

where ax is the path difference interval between samples. Using the transform pair 
relationship. 



(3) 


it is seen that sampling in the path difference domain is equivalent in the wavenumber 
domain to convolution of the spectrum by a series of delta functions spaced apart in 
wavenumber by multiples of ^ . 

Figures Ic and 2c illustrate (2) and 


i- 

P(a) * AX 



(4) 


respectively. 
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Sampling Theorem ORIGfWflL PABp ,e 

POOR QUAUn 

The sampling theorem has a deep significance in communication theory and other 
scientific disciplines involving measurements of signals at discrete intervals. It 
states that a bandlimited signal which has no wavenumber components above 
is uniquely determined by its values at uniform intervals less than 2aniax ^part. To 
capture the spectrum in Figure 2a that created the interferogram in Figure la, the 
sampling illustrated in Figure lb must be done at path difference intervals less than or 

equal to — cm. That is, 

‘^°^max 


AX < 


1 

^o^max 


(5) 


For the case of equality, the sampling rate is at the Nyquist interval. In (4) then, the 
spectrum is replicated at intervals of 2aff,gx» 


Truncation 


Of course, the sampling process cannot be continued for an indefinite period of 
time. If for no other reason, then's will be an upper limit on the number of samples that 
the system computer can process. Tnerefore, the series of sampling triggers is truncated 
by either the rectangular pulse shown in Figure Id or by^some other functional form that 
eliminates sampling for path differences in excess of t ^ from ZPD, Xt is the total 
path difference allowed by this truncation process. If the pulse in Figure Id is labelled 
A(x), then its transform is given by 


sin TraXt 


A(o) ■= Xj “ Xt S'nc(c'Xt) 


( 6 ) 


The sine function has its first zero at a = 


- 1 


X*~" 

the resolution Aa of the procedure. Figure 2§ depicts the sin x/x type behavior and in 
Figure 2e the convolution of the preceding signal with A(a) is shown. Normally, 
instead of th^. rectangular pulse, a tapered function known as the apodization function 
is used. A variety exist with each having attributes such as low sidelobes or fast 
rolloff that are attractive for particular applications. The particular form used for 
UVIS was described in PGW. For this study, the merits are not important; it is suffi- 
cient to use the rectangular truncation in the rest of the analysis. In summary, at this 
point in the development, the signal in the path difference domain can be described by 


and oftentimes this wavenumber is termed 
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a 


Figure 1. Fourier transform spectrometry signals in the path difference, x, 

domain, (a) The interferogram. (b) An infinite sequence of sampling 
triggers spaced ax path difference units apart, (c) The interferogram 
after sampling, (d) The apodization function (here represented by 
simple truncation), (e) The sampled interferogram after truncation. 





a 


Figure 2. Corresponding signals in the wavenumber, a, domain, (a) The wavenumber 
spectrum of the incoming light which is responsible for the interferogram 
produced by the spectrometer, (b) The Fourier transform of the sampling 
trigger sequence. In this domain, the delta functions are separated by 
the reciprocal of the sampling interval in wavenumber units, (c) The con- 
volution of the spectrum with the sampling function, (d) The Fourier trans- 
form of the apodization function, (e) The convolution of (a), (b), and (d). 
(f) The sequence of sampling triggers applied by the discrete Fourier trans- 
form to the output, (g) The sampled output spectrum which has a resolution 
given by the reciprocal of the total path difference displacement. 
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(7) 


where N 


y(x) “ y(x)a(x)A{x) » 

= • 

Ax 
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Fourier Transformation 

In Appendices A, B and Cj, there are appropriate derivations leading to a formulation 
of the discrete Fourier transform. When (7) Is operated upon by the transform, the 
discrete spectrum 


P(nAa) 




y(i<Qx;o“'^2Trnk/N 

N 

2 


n = 0,i,2...,N-l 


( 8 ) 


Is produced. Consisting of N samples (like the Input interferogram function), the spe 
tral resolution Aa is 


“ Nax (8) 

Therefore, the total range covered by the discrete spectrum Is 

ot “ ^Acr “ AX = 2amax 

If sampling is done at the Nyquist interval. However, as reported in Appendix C, this 
band contains the wavenumber samples 


N 

P(nAcr) for 0 < n < g 
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and the samples 


N 

P£(N-n)Ao] » P(-nAo) for 2 < n < N 

The actual spectrum is of course real only for positive wavenumbers. The spectrum 
recovered by the transform in the second ha'(f of the wavenumber span is the image of the 
spectrum ^nto negative wavenumbers because of the convoluticn shown in Figures 2a-c. It 
is a copy of the real spectrum and does not contain any information. All of the informa- 
tion of importance is contained in the first ^ samples produced by the transform. 

OTHER CONSIDERATIONS 


In this section several complications must be discussed that impact the operating 
characteristics of the spectrometer. These include the techniques employed to produce 
samples evenly spaced in path difference, the errors resulting from imperfections in this 
technique, and a discussion of the aliasing phenomenon. 


Path Difference Sampling Technique 

For UVIS, the interference pattern produced by the interferometer is sampled at 
fixed intervals of time using triggers produced by an external pulse generator. To get 
measurements of path difference, x, the spectrometer is also used to produce the inter- 
ference pattern for an incoming He-Ne laser beam. Being a coherent signal, the inter- 
ference pattern is a sinusoid that does not decay in amplitude as path difference is 
increased away from ZPD. And, since it is a monochromatic light source, the wavelength 
of its interference pattern fMist be equal to its electromagnetic wavelength. This sinu- 
soidal signal is then sampled at the same rate using the external pulse triggers. This 
rate must be sufficient to enable the accurate calculation of the location of some fea- 
ture of the sinusoid that can be used to monitor path difference. For UVIS, the zero- 
crossings were selected as that feature. After subtracting any bias present, the samples 
are used to compute the zero-crossings by means of a linear interpolation algorithm. If 
this could be done perfectly then ao each zero-crossing was reached during the production 
of an interferogram, it would be known that .3164 pm of path difference had been tra- 
versed. To determine the zero-crossing location to high accuracy requires a high number 
of pulse generator triggers per intirferogram fringe. Both the Interval between triggers. 
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PD, anu the period of the fringe pattern are based on time and can be related. The 
movement of the frrnge pattern wnh time depends only on the axial deformation of the 
heated aluminum bar as described In PGM. If a voltage is applied to the coll surrounding 
the bar such that Its length changes by HR ym/msec, then the number of msec per zero- 
crossing is given by . 3164/HR, Subsequent division of this by PD msec/samples yields 
the number of samples per zero crossing, R, Mathematically, 

R - ,3164/(PD . HR) 

PD and HR must be chosen so that sufficient samples are collected per zero-crossing to 
ensure that the zero-crossing locations are determined with minimum error. 

To explore this further, a sinusoid was sampled by use of a computer program for a 
range of values of R, Figures 3, 4, 5, 6 and 7 show the error between the Interpolated 
locations of 125 consecutive zero-crossings and the actual locations for R equal to 
approximately 6, 3, 2, 1,5, and 1,2 samples per zero-crossing, respectively. The maximum 
error in pm can be selected to characterize the effects of the sampling process. In 





Figure 3. The error in microns resulting from locating the zero-crossings in 
the laser reference interferogram with a linear interpolator and a 
sampling rate of six equal time spacing samples per zero-crossing. 
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Figure 8, this parameter Is plotted versus R. Also, the graph contains the ratio of the 
maximum error to the zero-crossing separation distance plotted versus R. For example, 
at a sampling rate of 1,5 samples per zero-crossing, the maximum error Is about 
1,85, 10“2 pm, which is about 5 , 9 % of the zero-crossing spacing. 

As discussed in P6W, the zero-crossings are used to determine the approximate values 
of the Interference pattern produced by the solar signal. These constant path difference 
samples would be unaffected by low frequency changes In the bar heating rate and also 
would be in the appropriate x domain for transformation by the discrete Fourier transform 
algorithm. The true cost of the location error discussed above cannot be determined 
without using these mislocated zero-crossings to trigger the sampling of the solar signal 
Interferogram. To simulate this, the following Fourier transformation relationship can 
be derived and conveniently used. If 

P(a) » e“n^(c"Co)^T2 . T (12) 

where ao 1s the wavenumber of this Gaussian function’s peak value, then 
1 

y(x) » 2 e cos 2TraoX (13) 


with 


T " l/7r(a'-ao) 

0* 1s defined to be the wavenumber at which the Gaussian function's amplitude is 1/e of 
Its peak value. The path differences, x, are In microns. If (13) is sampled at the 
zero-crossing points determined from a linear Interpolator and which therefore are 
mislocated, then a set of samples Is produced which Is useful for this evaluation. When 
used as the input data to the Fourier transform algorithm, the function shown in Figure 9 
was produced. The solid line represents the original spectrum given by (12) with oq “ 
11,500 cm“^ and o' = 11,000 cm“^. The errors introduced Into this sampling and transfor- 
mation procedure are those illustrated in Figure 6. This case, which has a maximum 
mislocation error of 1,85 • 10**2 results in a recovered spectrum that is most 
marked by the introduction of a non-negligible artifact at about 4,350 cm“^. Wi^-hin 
the skirts of the recovered Gaussian, there is a small shift towards higher wavenumbers 
as well; The impact of these spectral errors will be discussed later in this section. 


n 
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Figure 8. The variation of the maximum xero-crossing mislocation error in 
microns with the time interval between sampling triggers, PD, in 
msec. The additional axes show the variation of this error as a 
percentage of the laser reference zero-crossing path difference 
.3164 pm as a function of the number of samples per zero-crossing, R. 
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NYQUIST LIMIT 


MAGNITUDE 



Figure 9. A normalized spectrum produced by the FFT with and without 

( } the interferogram samples being mislocated by the errors 

plotted in Figure 6, 


Interferogram Centering 

The only mechanism that has been found to be capable of smoothly moving the inter- 
ferometer optics by the small path difference required for visible and near-ultraviolet 
spectral measurements is the heated rod. One might suspect that this technique would be 
flawed because of inordinate amounts of time that would be required for the bar to be 
heated and cooled through the appropriate temperature range. With aluminum as the 
selected material, this is not the case. For a reasonable voltage, the Nichrome heating 
coil can easily create a bar expansion rate for heating of around 10 pm/sec. A cooling 
rate nearly as fast can be produced by using an auxiliary fan. If, for example, a spec- 
tral resolution of 100 cm“^ were desired, then from (9), equals 100 pm. For the 
above expansion rate, a complete scan can be obtained in about 10 sec. For this length 
of time, the atmosphere can certainly be considered steady-state. 
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However, the true cost of the heated bar technique Is the lack of registration. 
Consecutive scans cannot be started at the same path difference. By use of mechanical 
stops, they can be started within a few microns of the same path difference. Therefore, 
it Is necessary to collect more samples triggered from the pulse generator than are 
actually needed for the interfercgram. By starting a little in advance of the desired 
scan starting point and by finishing a little beyond the desired end point, there will be 
a cushion of samples that can be discarded as needed to bring an interferogram into proper 
registration. If M is the number of pulse generator triggers and if aM is the additional 
samples in the cushion, then the number of zero crossings is given by (M+aM)/R. The 
situation can be further complicated if a path difference sampling interval ax different 
from the zero-crossing interval is desired. This can be done in the software by linear 
interpolation between zero-crossings or by skipping integer numbers of zero-crossings. 
If the number of final, or constant path difference samples per zero-crossing is K, then 
the total number of interferogram samples is 

1 

(N+aN) = (M+aM) • r • K 04) 

where aN is the additional numbers of constant path difference samples that is produced 
by the aM samples in the cushion. Zero path difference is found by visual inspection 
of the final interferogram. Then the ends are truncated as required to produce the proper 
Interferogram length. 

Of course, the interpolation and the imprecise determination of zero path difference 
give rise to potentially important phase errors. It is for this reason that the two- 
sided interferogram and discrete Fourier transform are used in UVIS. By taking the 
magnitude of the computed spectrum, these phase errors can be avoided. Consider the 
following explanation. Pe{a) is the actual spectrum which is even (i.e., symmetrical 
about ZPD), ,j)(cr) is the function describing the phase errors present. Then the 
observed spectrum is 

= p(a)-iq(a) = J dx (^^) 

.*00 

By taking the magnitude using 

the actual spectrum is recovered regardless of the form of (j)(cr). 
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Aliasing 


The final additional consideration that is necessary for describing the operating 
characteristics of UVIS is aliasing. This phenomenon will be described using Figures 10 
and 11, which are modifications of Figures 1 and H. In Figure 11a, the spectrum from 
Figure 2a is shown after it has been attenuated on the low wavenumber side in a fashion 
that is typical of the effects of an optical filter. If its total bandwidth is now less 
than 1/2 of the upper wavenumber cutoff, amax» benefits of aliasing can be 
reaped. Figure 10b shows constant path difference sampling performed at twice the 
interval used in Figure lb. This results in a replication of the spectrum at intervals 
of amax» snown in Figure 11c. Then, = Nax is twice the total path difference 
as before, and the spectral resolution is halved as shown in Figure Ilf. For the same 
number of samples, N, the computed spectrum can be found in the first N/2 spectral samples 
as shown in Figure llg. The example shown produces a backwards replica of the spectrum 
because of the wavenumber domain convolution process but this is easily corrected by the 
software package for UVIS. 

The governing equation for aliasing is a modification of (5). If ax U, selected 
such that 


1 

1 2{<7max-«^o) 


(17) 


where oq is the lower wavenumber cutoff of the filtered spectrum and where 


f^max “ <^o " c^max/'^ 

then aliasing can be a beneficial attribute of the spectral recovery scheme, m is 
an integer and is called the spectral order. For the example above, m = 2 and ao 

2. cymax/2. 

One specific requirement for this technique to work is that there are no spectral 
components below oq* there are, then the aliasing will introduce these components 
as contaminants to the replicated spectrum produced by the convolution. Therefore, the 
artifact produced by the laser beam interferogram Interpolation procedure that is shown 
in Figure 9 could not be tolerated if aliasing was being used. 
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Figure 11. Corresponding signals in the wavenumber, a, domain, (a) The spectrum of 

the incoming light is filtered so that the low wavenumber portion is missing, 
(b) With &x twice as large as in Figures 1 and 2, the sampling function 
triggers in the wavenumber domain are spaced half as far apart, (c) The 
convolution of (a) and (b). (d) The Fourier transform of an apodization 

function having twice the total path difference displacement as before, (e) 
The convolution of (a), (b), and (d). (f) The sequence of sampling triggers 

applied by the FFT to the output, (g) The sampled output spectrum with 
twice the resolution as in Figure 2 but with the same number of samples 
because of the use of aliasing. 
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From all of the preceding discussions, a system of equations has been developed 
which completely describe how an instrument system such as UVIS will operate. These are 
recapitulated below for review. 


1 

2(amax'«^o) 

(19a) 

"’(omax“cfo) * cJmax 

(19b) 

i 

N + aN = (M + aM) . R . K 

(19c) 

N = 2^, J an integer 

(19d) 

Tt = (M + aM) . PD 

(19e) 

Xt = (N + aN)ax 

(19f) 

1 

Ao = Nax 

(19g) 

<?t “ ^'Omax 

(19h) 

R = .3164/{PD . HR) 

(191) 


These have been recombined in a simple computer model to produce values of omax* ax, 
K, N, M + aM, Xt, and Tt given input values for amax* omin» Ao* PD, HR, and aN. 

An initial value for amax is entered and the output amax value produced is the wave- 
number needed to satisfy (19a,b,d,g) in combination. The Fortran code for this model is 
included as Appendix D and is suitable for use with the DEC RT-U version 3b operating 
system. 
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MODEL RESULTS 


The following results have been produced using the model subject to these conditions. 
First, the cushion aM has been set at 10 % of N. Secondly, the bar expansion rate of 
10 pm/sec is fixed. However, the pulse generator triggering rate PD is variable. 
Based upon the result given in Figure 9, it seems prudent to require that PD be less 
than 20 rnsec/sample. This is equivalent to saying that R should be kept greater than 
1.5 reference samples per zero-crossing. 

One of the input parameters is ao» the output spectral resolution. Of course, it 
would be desirable to keep this as small as possible. In the visible and near-ultraviolet 
wavelengths, most scientists and engineers tend to use wavelength, ^, as the unit of 
measurement rather than wavenumber, o. Therefore, it is Important to point out that a 
constant resolution Aa does not correspond to a constant resolution a\» The wavelength 
at which the resolution element Aa is located influences the corresponding resolution ax. 
This is illustrated in Figure 12. For example, a resolution of 100 cm”^ at 1.0 pm is 
equivalent to a wavelength resolution of lOOA; but at 0.3 pm, it is equivalent to 

9A. Aa = 40 cm“^ is the same as 40A at 1.0 pm and 3.6A at 0.3 pm. As will be 

seen in the model results, quite high wavelength resolution in the ultraviolet is possible 
if sufficient samples can be collected and transformed. As seen in (19g), ao is dependent 
upon AX, the interval of constant path difference. It should be observed that the 
original sampled data contains higher resolution information because of the requirement 
that there be enough time samples per zero-crossing to make possible accurate interpola- 
tions of the zero-crossing locations. This information is lost, however, because of the 
FFT algorithm requirement for evenly spaced data points in path differences, 

Before specifying the optimum operational characteristics of the UVIS, it is 

instructive to exercise this model to examine the variability of the important output 

parameters for a range of input parameter values. Figure 13 shows how M and N vary for a 
number of amax values. Only spectral order 1 (m=l) solutions are plotted and the spec- 
tral resolution is 40 cm“^. The pulse generator trigger period, PD, is 20 msec. The 
discontinuities in the curves are the result of the requirement that N is a power of 2. 
Within each plateau region, the value of K, the number of output samples per zero-crossing, 
is also given. 

In Figure 14, the effects of varying PD is examined for the different spectral 
resolutions, 100 cm“^ and 40 cm"l. Again only solutions for spectral order 1 are shown 
and is set at 30,000 cm”^; K is constant for all solutions and equal to 2 samples 
per zero-crossing. 
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max 


(cm - ) • 10 


The variation of M and N with maximum wavenumber. Only spectral 
order 1 solutions at a resolution of 40 cm“'^ are shown. PD is 
20 msec and the number of samples per zero-crossing, K, are shown 
in-between the discontinuities. 





UVIS OPERATING CONDITIONS 


To speclftcally identify the operating parameter values for UVIS, it is necessary 
to identify the appropriate values of oo» f>otential applications of 
UVIS. As discussed in PGW, the instrument uses silicon photodiode detectors that are 
responsive from about .25 ^m - 1.1 jjm. That is, they are sensitive to radiation from 
the near ultraviolet through the visible and into the near-infrared. Figure 15 is the 
solar spectrum plotted against wavenumber as reported by Mecheri Runnel (1980). Tenta- 
tively, this can be divided into three broadband regions that can be used for the follow- 
ing applications. From about 10,000-15,000 cm"^, there are three water vapor absorption 
band systems. By monitoring the difference between the signal strength at the center of 
these bands to the signal level at unattenuated wavenumbers just outside these bands, 
estimates can be made of the amount of water vapor present in the atmosphere. Prelimi- 
nary results from UVIS for this application have been reported by Parsons (1983). From 
20,000 - 30,000 cm~l, the solar spectrum becomes increasingly attenuated by ozone as 
the wavenumber increases. By monitoring the gradient of the sun's intensity with wave- 
number and hopefully by monitoring its value at specific ozone absorption lines, UVIS 
may be able to measure total ozone amounts. Finally, the spectral band between 15,000 
and 25,000 cm"! encompasses the visible wavelengths and an interferometric spectrometer 
operating within this band would obviously have many uses. 

Therefore, these three bands are adopted for this modelling effort to give reason- 
able values for amax ^^min* Aa cannot be specified at the desired levels because 
of computer memory constraints. The LSI-11/02 system used in UVIS has a memory capacity 
of 32K words of which 28K words are available to the programmer. This limits the value 
of N to 512 or less and M to ISOO or less, ao will be selected at the highest resolution 
possible to accommodate this constraint. Figure 14 .shows that at ao “ 100 cm"^, for 
example, the required number of original samples, M, rises dramatically for smaller 
values of PD to produce the same number of final samples, N. 

To stay within the memory limit and also to keep PD less than 20 msec to avoid the 
Interpolation artifact problem may necessitate a very narrow range of operating condi- 
tions. This is examined in the following illustration. Figure 16 shows the variation 
of N and M with ao for the bandpass covering 10,000 - 15,000 cm"^ and for R = 1.5. This 
graph shows that for this band, both spectral orders 1 and 2 are usable. Of course, m=l 
Implies that the resultant spectrum covers the entire range of wavenumbers from 0 cm”^ 
to 15,000 cm^^. Actually a ,y,ax 15,802 cm“^. For m=2, oq = 7,901 cm“^ instead of 
0 cm"l. To conform to the computer memory limit, the maximum resolution that can be 
achieved for R = 1.5 is 61.73 cm”l. Table I summarizes the information contained in 
Figure 16 and adds similar Information for the same bandpass but for the R = 2.0 and 
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Figure 15. The solar radiance spectrum from Merrharikunnel and Richmond (1980). 
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AO" (cm ) 


The variation of M and N with spectral resolution for a spectral 
bandpass of 10,000-15,000 cin“l and 1.5 samples per zero-crossing 

, N (spectral order 1, K=l)j , 


M; 


(spectral order 2, K=l/2). 







R * 3.0 cases. For R * 3.0, it indicates that 123.46 cnr^ Is the maximuni resolution 
achievable. 

Table II applies to the ultraviolet bandpuss containing the ozone absorption fea- 
tures, For sampling rates of 1.5 and 2,0 final samples per zero-crossing Ac = 61.73 
cm“l is possible but only for spectral order 2. For R = 3.0, again only 123,40 cm”l 
resolution is achievable. Finally, the visible wavelength bandpass is studied in Table 
III, For c,nax “ 25,000 cm*^ and oq h 15,000 cm"^, the possibility of using aliasing 
and spectral order 2 is gone. Unfortunately, for all three sampling rates the best 
resolution possible is 123.46 cm“^. 

To surnsnarize these findings Table IV is presented. It contains the pertinent infor- 
mation from these model results. When several operating modes are possible, the one 
with the best resolution and the highest sampling rate R is the one entered in the table. 
The Information entered for the 15,000 - 25,000 cm"^ case is also appropriate if the 
entire spectrum to 30,000 cm"^ is to be recovered. 


table I. SAMPLING RtQUiKEMENiS FOR WAVENUMBER REGION 800U-15000 cm’^ 


Ac Range 


N 



M 


(cm"l) 


m=l 

m=2 

R=l,5 

R=2.0 

R*3. 0 

15.43 - 

30.87 

2048 

1024 

3414 

4649 

6824 

30.87 - 

61.73 

1024 

512 

1704 

2276 

3414 

61.73 - 

123.46 

512 

256 

854 

1136 

1704 

123.46 » 

246.92 

256 

128 

426 

569 

854 


TABLE II. SAMPLING REQUIREMENTS 

FOR WAVENUMBER REGION 20000 

-30000 cm 

Ac Range 

N 



M 


(cm"l) 

m=l 

m=2 

R=1.5 

R=2.0 

R=3.0 

15,43 - 30.87 

4096 

2048 

3414 

4549 

6824 

30.87 - 61.73 

2048 

1024 

1706 

2276 

3414 

61.73 - 123.46 

1024 

512 

854 

1136 

1704 

123.46 - 246.92 

512 

256 

426 

569 

854 
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TABLE III. SAMPLING REQUIREMENTS FOR WAVENUMBER REGION 15000-25000 cm"l 


Acr Range 

N 



M 


(cm“^) 

m=l 

m=2 

R*1.5 

R=2.0 

R=3.0 

15.43 - 30.87 

4096 

CO 

4551 

4551 

6827 

30.87 - 61.73 

2048 

1706 

2274 

2274 

3412 

61.73 - 123.46 

1024 

854 

1138 

1138 

1707 

123.46 - 246.92 

512 

426 

568 

568 

852 


TABLE IV. SUMMARY OF OPTIMUM OPERATING CHARACTERISTICS IN ALL 
THREE BANDS. 


N 
M 

Aa 
R 
m 

AX 

Xt 
Tt 

%ax 
PD 

CONCLUSIONS 


a band 

10-15,00 cfH-1 

512 

1137 

61.73 cm“l 
2.0 
1.0 

.3164 pm 
.3593 rnm 
17.033 sec 
15802.78 cm“^ 
0.0 cm"^ 

15.0 msec 


a band 

15-25,000 cin-1 

512 

851 

123.46 cm-1 

3.0 

1.0 

.1582 pm 
.1347 mm 
8.5167 sec 
31605.57 cm’^ 
0.0 cm"^ 

10.0 msec 


a band 

20-30,000 cm-1 

512 

1137 

61.73 cm~^ 

2.0 

2.0 

.3164 pm 
.3593 mm 
17.033 sec 
31605.56 cm“^ 
15802.78 cm"^ 

15.0 msec 


Using the design of the ultraviolet interferometric spectrometer (UVIS) and the 
fundamentals of the discrete Fourier transform, the optimum operating characteristics of 
the instrument have been developed for three broadband spectral regions. Each of these 
regions contains the signatures of atmospheric constituents which are important to several 
appl ications. 
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APPENDIX A 


ORIGINAL PAGE IS 
OF. POOR QUALITY 


THE FOURIER TRANSFORM OF A GENERAL PERIODIC FUNCTION f(x) 


From Appendix B, it is seen that 
F {1} = 6(a) 

One of the useful general properties of the Fourier Transform is frequency shifting or 

I 

modulation. Mathematically, this is evident as the following derivation illustrates. 

CP 

F {e^'^TiooX f(x)} = f f(x)e‘i2TT(a“ao)'< dx 

JucO 

= f*(c-ao) 


where f is the Fourier transform of f. If f(x) Is the constant 1, then 


F = 5(a-ao) 


This can be usefully applied to the following expression 


f(x) = 


n=-oo 


i2TrCt") X 


which states that any periodic function can be represented as the infinite sum of comp! 
sinusoids weighted by coefficients given by 


i_ f'"" - 

Xq J ^(x)e 

^-Xo/2 


-i2Tr(j, )x 
^0 


dx 


Using the Fourier transform definition. 


-12-fraX 


F{f(x)) = 


/ 


M 


12n(x )x 


y ane 


n=-co 


dx 


B 


00 



ns-o 



-i 2 tioX 

e dx 


Using the frequency shift property, 


F{f(x)j = 



6 (o 
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APPENDIX B 


ORIGINAL PAGE fS 
OP POOR QUALITY 


THE FOURIER TRANSFORM OF A CONSTANT 


Using e"nl^i as a convergence factor, 

F{1} = lifii j e-nlxl dx 

n-)-o 


= 1 im 
rr»-0 L 

f e(n-’ 12 -rr 0 )x dx + 

-00 

f 

■'0 

e-(n+i2iT0)x dx 

= lim 
n+0 

exp(n-i2TT0)x ° 
n~'i27rcr _oc 

) 

00 

exp(-(n+i 2 iT 0 )x) 

n+i2Tr0 ^ 


r 1 

1 1 



= 1 im 
rr^O 

n”i2ir0 n+i2iT0 




n+i2TT0 + rr*i2ir0 1 


lim 

n-K) 

2n 

= 1 im 
rr»-0 

n^+(2iT0)^ J 

= 



In the limit, there is a singularity at 0 = 0, By invoking L'HospitaTs Rule, the limit 
can be evaluated, however, and found to be 


F{1} = lim ( 2 r,) 

n+0 


0 


at 0 = 0 
at 0 0 


The area under the function can be obtained from 




mCO 


2r) 

n^+(2n0)^ 


da 


2n dw 

** /.cx, 2n 
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where w = Snc has been used as a substitution. Proceeding, 



2n 



n2(H^^2) 


dw 

2ti 


00 




ndx 

1+X2 


using x = - as a substitution. Then, 
n 


00 

i I dx 

IT i+x2 


i 

TT tan"^ X 


oo 


l«»00 


j, 

ir (tan”'*' “»-tfln‘*^(“<») ) 


JL JL JI2L 
“ n (2 - (2 )) = 1 


In the limit then, the Fourier transform of a constant of magnitude 1 has been shown to 
have infinite amplitude and unity area at wavenumber a = 0 and to be zero elsewhere. 
This is the definition of the delta function. Thus, 


F{1} = 6(a) 
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APPENDIX C original PAGIii ’ m 
OF POOR QUALITY 

DERIVATION OF DISCRETE FOURIER TRANSFORM 


Starting with equation (7), 


y(x)a{x)A(x) = 



= y(x) 


there is the need to compute its Fourier transform. Because the data is known at dis- 
crete points and because of the computer's obvious advantages for doing complex arith- 
metic manipulations, it is advantageous to compute the transform at discrete points. 
The discrete Fourier transform equation is derived as follows. If ' y(x) '•s the x 
domain function to be transformed, we can use the expression in Appendix A. The an 
coefficients can be calculated as 




g)-' 

E Y(x)6(x-kAX) 



n 

-i2Ti(^ )x 
c 
e 


dx 



So 



ORIGiNAL PACr: 
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F{y(x)} 


r-n 

L 


6(a 




n is the index of variation In the spectral dcir.ain and the delta functions are separated 
by wavenumber intervals of ^ 

Discrete Fourier transforms, when implemented in digital computers, can be made much 
more efficient if particular parameter choices are made (Chamberlain, 1979, pp. 315-325). 
First, if Xc is set equal to X^, the total path difference traversed by the inter- 
ferometer, then 


M M 1 

*0 ■ ° N 

where N is the total number of interferogram samples. Then 


F{y(x)} = 



n=-oa 




The an coefficients are sampled at intervals of wavenumbers. The discrete Fourier 
transfOMti is addit’unally simplified if the range that n is allowed to assume is limited 
to N values. The impact of this decision is analyzed in the following calculation. Let 
n = r. 


C2 
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“ P(r6a) 




id. 

-iZrrrk/N 

Y(l<Ax)e 


6(a - X*^) 


If n « r + N, 


1 

>'{Y(x)}n=p+N = P((r+N)Aa) = 



Jk 

“127r(r+N)N 


(r-fN) 
s(o - ) 


But 


e-12,(r+N)k/N =, e-i2;.rk/N ^-IZuk = e-i2irrk/N (<.os 2,,k-isln27tk) 
= e-i2j[rk/N fgp |( integer 


F{7(x))„=rH,N 


X, 


(?)■' 

E-i2irrk/N 

Y(kAx)e 6 (a 




k=- ^ 
K 2 


That is, at a = (r+N)/\cr, the discrete Fourier transform yields the same result as for 
a - TAcj. If N samples are examined, all of the information available from the 
transform will be recovered. 

The discrete transform equation then is given by 



n = 0,1,2,, ..,N-1 


P(nA0) 


Y(kAx)e"i27rnk/N 



This replication also holds for n » r-N and for positive and negative multiples of N. 
That is, * 


P(riAo) • PC(n + JN)Aa3 d « 0, ± 1 , ± 2 ,... 
In particular, it is important to note that 
P[(N-n)Aa3 “ P(-nAa) 

That is, the computed spectrum from wavenumber samples 
lies between and 0, 


N 


I 


to N is the same spectrum that 


ORIGINAL PAGE 511 
OF POOR QUAUTV 


C4 



ORIGINAL PAC'!i B 
APPENDIX D OF. POOR QUALITY 

SIMULATION PROGRAM 


C 

- r. 

PROGRAM FTSSIM/fR 

yvRiT.TEN BY C> L. ^.PARSONS ( CDDE_.9 7 1) ON 10 / 1 «/H2 

c 

r 

DESIGNED^ TO 

COMPUTE PERFORMANCE CHARACTERTSTTCS OF A 

C 

- C 

FLEXIBLEf COMPUTER-CONTROLLED FOURIER “TRANSFORM SPECTROMETER ' 
GJ..V.EJil_Tri£ FOLLOWING INPUTS!. .. . ... . ' 

n 

SIGMAX 

- UPPER WAVENUMBER CUTOFF OF INCOMING LIGHT 
SIGNAL SPECTRUM .(IN, ,CM-1) 

c 

CL 

SIGMIN 

DELSIG 

- LOWER wavenumber CUTOFF (IN CM-l) 

- DESIRED SPECTRAL RESOLUTION (IN CM-11 „ _ 

C 

r 

PO 

- EXTERNALLY CONTROLLED INTERFEROGRAM 
SAMPL1NG..PERI0D (IN. MSEC 1 

c 

r 

HR 

- HEATING RATE OF OPTICAL TRANSLATION 
me c H a N I SM_LIN Z.ER0.«.Cii.a3 SliJGSZMS E CJ 

c 

-,.C. 

XBAR 

- LASER REFERENCE SIGNAL ZERO-CROSSING PATH 
OJFFERENCE = -316A MICRONS 

c 

- - .r 

DELN 

- NUMBER OF FINAL SAMPLES NEEDED TO CUSHION 
... T.HE...CENJFRTNR OF INTFRFFRnRRAM (TN PERCENT 

c 

r 


OF N) 

c 

THE OUTPUT PRODUCTS ARE! 

(1 

_...3I.GMAX1 

- UPPER. WAVENUMBER CUTOPP OF MEASURED 

c 


SPECTRUM (IN CM-l) 

r 

- M 

......SPECTRAL ORDER 

C 

OELX 

- INTERFEROGRAM SAMPLING INTERVAL (IN MICRONS) ' 

(' 

YFXTRA 

- JlQa_D.Fl_F.INALXS.AMP,L£SZjJERQ=.CPDSSING 

C 

N 

- MO, OF FINAL SAMPLES USED IN FFT PROCESS 

. r. 

- MTOTAL 

- MO- OF ORIGINAL SAMPLES NE.EDECLJO PiLODUCE N 

C 

XTOTAL 

- TOTAL TRANSLATION DISTANCE NEEDED (IN MM) 

r 

- ._ . TTQTAL 

.J OTAl. TIME REtJUIREl) FO R TRAN S L A.1 1 0 fiULI N SEC) 

C 

R 

- NO. OF ORIGINAL SAMPLES/FINAL SAMPLE 

.£ 

RPRIME. 

-..NO, OF original SAMPLES/ZERO-CROSSING 

C 

_Ji£AL_Ml.tlIAL 



C' 


■ 10 TYPE " ENTER TH E A P PR p.K I.M ATE UPPE R WAVENUM B t R_C UJ QFTLJ 

READ FREE (UrERRslO) SIGMAX 

— L5 _JDlE£__L-£!yJ£.R-.IH E. . A P.e R.UX LW A.I E . L 0 w , ^AV £N tj m EELi: UXO,F_E__!. 

REAP FREE (11,ERR = 15) SIGMN 

^aCL-l-Y,g.£,JL-£XaEfX-I.ri£. D £S I RE D,..XX1IAL„0JJ T £UJ_5PXCIML.RES.0LUXL0.N__C£ d- 1 
READ FREE (11»ERR=20) DELSIG 

L.Y fiE_.X..tiiE_£ALCU LA TXaN S .. A S S.U mE_J H AX. A_.h£^ViL_L A S E R_13_U5 ED 

TYPE ' THIS implies A REFERENCE ZLRU-CRUSS ING PATH DIFFERENCE * 

-TYPE * OF MiC PUNS...! 

i"RIIE (10,25) 

-J>5....E:uR>iiA.T...-L LLL I T...I S ,.A LS 0 .A5S.UMED._m AJ..EXXRA_BAMpXL3_.ARJ.LNf / ) ' 

TYPE * TO ENSURE THAT A CENTERED HLOCK OF DATA CAN BE EXTRACTED ♦ 


YPE_!. JiKQ ‘i_t HE_K A vv.. DATA 

TYPE ' THE default VALUE 


,e.ILE_X.._ 
FOR THE 


NUMBER OF EXTRA SAMPLES IS » 


r.tP.L_L_LO„P_ERCENj_uF.. THE J;.INaL. SAi^PLES NEtl^FD AS CUMPUIED 6EL0N * 




0 RK 5 INAI PAGE FG 


XTOTALsMrorAL*OELX)»lF -3 


OF POOR QUALITY 


^IIi 1 IAL!^ 1 IQXA 1 .AP.O>U£.* 3 -,.= -, ,,, 

rt 3 MrUTAL/((l.»+(DfcLN/l 00 .U) JAN) 

«P.aiiiE:saftX£XIRA-„ _ 

c 

c- £RXxa~QaT..,R ESULia^cLjaaLxiL t 

c 




LO.lL^.QHiiA.'UXlIlD 

105 FORMAT ( ' FOR A DESIReO WAVENUMBER WINDOW FROM %F 10 , 2 f 

alUILM - LJCILJLiXlfl.. a *I^C M- L.JL ,) 

loa FORMAT ( » (OR SF8,2#'A TCJ SF8,2,M). ' ) 

UJ..EORMaT_J_»_A.NO_A^.SPECT,RAL.RESOLUTIUN of »,F 7 . 2 r' CM-l » ) 

112 FORMAT ( ♦ (UR %F 7 . 3 #'A AT SOOOA) » ) 


...,ajJ5..F:LmMM-,..LJL,XHE.XaLL0wIN ,F0URiE.R.JRA£i5EiJR.M-.SRECXROI:iEriR.„DE8lGiv.S ), 
120 FORMAT ( * ARE R08SIBLE: *///} 

,.-iafa-.FQHM..A.C-.C„-l^(.rHE LASER fiEFERENClE_.C.HANNj:L„2ERO-rCRjO.SS 1NG*J’AXH.J.«) 

130 format ( * DIFFERENCE IS ASSUMED TO BE .31b^ MICRONS ) V/J 

.....iJlS ^E01iMAL..(_!^AKD W A RE - V AR lA B LES JZJ 

lao format ( ' THE PULSE GENERATOR SAMPLING PERIOD = SF5.W' MSEC *) 
^_l<l5-.TJjHi^.AT„ (— 1-J.OE JAR SEATING. RATE.- = .J.*XA.3fX..Z£R0--CR0S$Ii'J6/MSEC 
150 format ( » SOFTWARE VARIABLES — V) 

— LS5-£QiiMiU_(^X-^0E„,SPECJRAL OROER.J./.II 

160 format ( ' THE OUTPUT SPECTRUM UPPER CUTOFF = %F10*2, 

*• rM»i » > X j 

162 format ( ' the output SPECTRUM LOWER CUT0FF”%F 1 0 .2, » CM-l ' J 

— L(tx5_Fil&aA.L_X_: I.u£-JlAIH..DXEFiLHEj'ICE«INTi:aVjiL-.5™X<L£Z^^.,JL«M.lCRJ)RS_L^ 

170 FORMAT ( » THIS REQUIRES A SAMPLING RATE/ZERO-CROSSING = 


-AF-6 

175 FORMAT ( ' THE NUMBER OF POINTS NEEDED FOR THE FFT = M5/J 

X8.0_4?aftMAX^-j: A.SJUMXAlG-A--CUSHIUNLjllUAL-.tQ_JL*ES,J_» 

n» PERCENT OF THIS NUMBER OF POINTS^ ' ) 

XaS-EOlifcULL-C— A-LQ t A L-DF *.1. JJiil GIW AL_ JAMBL E8_LSL_N E£ D ED .JJ. ) 

190 FORMAT ( ' THE TOTAL TRANSLATION PATH DIFFERENCE = %F7,a, 




195 FORMAT X • AND THE TOTAL TIME NEEDED/SCAN = SF8.4»' SEC*) 

afliL.. F .QJM AI ^. ( ^.. * DIE. JsllitiaEJL-flE-.Qfi.IG II!iAL„.aAMPL£aZ£ IN AL_S AMP L£ * ^ 

*F8,A) 


2.ttS-JilfiMAI_X-J: AKDLj:B£-NUm£JLJ3£-.QLSXG.I-N A L S AM PL.£Sy.Z£fiflr.CJCI.aSJ.NQ-.£Xt 

*F 8 ,‘l) 

^R.IXE, ..( 1 2^.10JU 

WRITE (12#105) SIGMIN,3IGMAX 

viUiWsuj^cuaxGMX N 

wLMAXs(U0/3IGMAX)*1£8 

-iviUX£-.(JL2^QA)«J<LMlN,.WLMAX _ 

WRITE (I2rll0) OELSIG 

UEUtLsflELaiGA ^5.—,, 

WRITE (12,112) OELWL 

,_j*j A11£-CX£UJ. 5 J 

WRITE (12,120) 

i'RXXE-Xl£,J.25J-. 

write (12,130) 

/iRj;.l£.,Xl£,.i35i 

WRITE (12,1'40) PO 

WR.LLL...,C 1 ?.U.A5.L.hR 

.'-RITE (12,150) 

i\RIXL. ( 12 « 1 55J-.M 


D 3 



ORIGINAL PAtiK S3 

wwne (i2#i60) siGMAxi OF POOR QUALny 

^RJLL£-U,?.<.162J^ 3 IG«JLW1 

WRITE (12#16S) OELX 

— uRU£-JJLa*l2jQJ JtEXIHA 

WHITE (12#175) N 

WRITE (12»I05) MT0T*aT 

flilI.T.£..-a 2 1.9.flLl,XiaiAi 

write (12,195) TTOTAL 

WRU.E^i2.fc£Q^LJt— 

WRITE (12,205) HPRIME ' 


L FOR M.GT.l LOOP BACK FOR OTHER ORDERS 

. ... C 

IF (M.GT.l) 60 TO 300 

GO TO AO O 

300 MsM-l " "" 

G.U -IQ.,.6.0 . „ 

C 

C a£ClD.£...IF QJ±LE£J£AS£.S_AR £-T_Q-HE RU N 

C 

— -_a0CL.XYPE_!..aF.aiTHEK.CASES , ARE_TO BE.-RUN,_T.YPE.JL, OTHERWISE. 0 / 
READ FREE (ll,EHRsA00) INUH8 

IF—UNLIM a^£D^(lT,_G a .TQ 9 ^0 ___ 

GO rO 10 


G WRAPUP 

-C 

800 TYPE * A SPECTRAL ORDER LESS THAN 1 WAS COMPUTED - ABORT ' 

.9hXLXaN.TXWUE 

STOP 

PJUD._ 








